Rational Positive Systems for Reaction Networks

Extended abstract The purpose of the lecture associated with this paper is to present problems, concepts, and theorems of control and system theory for a subclass of the rational positive systems of which examples have been published as models of biochemical cell reaction networks. The recent advances in knowledge for the genome of plants, animals, and humans now lead to increased interest in cell biology. Knowledge is needed on how a cell as a functional unit operates biochemically and how the reaction network is influenced by the genome via the enzymes. In principle it is possible to model the complete biochemical reaction network of a cell though this program has so far been carried out only for small compartments of such networks. Mathematical analysis for such reaction networks then leads to a system of ordinary differential equations or of partial differential equations. Often the ordinary differential equations are of polynomial or of rational form. The number of reactions in a cell can be as high as 15.000 (about half the number of estimated genomes) and the number of chemical compounds as high as 20.000. A detailed mathematical analysis of a mathematical model of the complete cell reaction network may therefore not be possible in the short run. Hence there is an interest to develop procedures to obtain from high-order mathematical models approximations in the form of low-order mathematical models. The formulation of approximate models requires understanding of the dynamics of the system, in particular of its algebraic and graph-theoretic structure and of its rate functions. It is the aim of the author to contribute to this research effort. In this lecture attention is restricted to mathematical models for biochemical cell reaction networks in the form of rational positive systems. These systems are called positive because the state vector represents masses or concentrations of chemical compounds and the external input vectors represent inputs into the network of externally available chemical compounds and of enzymes produced by the nucleus of the cell. The dynamics of the system is often modelled as a polynomial map but in this lecture attention it is restricted to rational maps (each component equals a quotient 2 Jan H. van Schuppen of two polynomials). Such a dynamics arises for example in the model of Michealis-Menten kinetics due to a singular perturbation of a bilinear system. The mathematical model of the glycolysis of Trypanosoma brucei is phrazed almost entirely …